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EXPERIMENTAL  STUDY  OF  COMBINED  FORCED  AND  FREE 
LAMINAR  CONVECTION  IN  A  VERTICAL  TUHEI 
By  Theodore  M.  Hallman 


SUMMARY 

An  apparatus  was  built  to  verify  an  analysis  of  combined  forced  and 
free  convection  in  a  vertical  tube  with  uniform  wall  heat  flux  and  to 
determine  the  limits  of  the  analysis.  The  test  section  was  electrically 
heated  by  resistance  heating  of  the  tube  wall  and  was  instrumented  with 
thermocouples  in  such  a  way  that  detailed  thermal  entrance  heat-transfer 
coefficients  could  be  obtained  for  both  upflow  and  downflow  and  any  asym¬ 
metry  in  wall  temperature  could  be  detected. 

The  experiments  showed  that  fully  developed  heat-transfer  results, 
predicted  by  a  previous  analysis,  v/ere  confirmed  over  the  range  of 
Rayleigh  numbers  investigated.  The  concept  of  "locally  fully  developed" 
heat  transfer  was  established.  This  concept  involves  the  assumption 
that  the  fully  developed  heat-transfer  analysis  can  be  applied  locally 
even  though  the  Rayleigh  number  is  varying  along  the  tube  because  of 
physical-property  variations  with  temperature. 

Thermal  entrance  region  data  were  obtained  for  pure  forced  convec¬ 
tion  and  for  combined  forced  and  free  convection.  The  analysis  of  lam¬ 
inar  pure  forced  convection  in  the  thermal  entrance  region  conducted  by 
Siegel,  Sparrow,  and  HalLman  was  experimentally  confirmed.  A  transition 
to  an  eddy  motion,  indicated  by  a  fluctuation  in  wall  temperature,  vras 
found  in  many  of  the  upflow  runs.  A  stability  correlation  was  found. 

The  fully  developed  Nusselt  nvunbers  in  downflow  were  below  those  for 
pure  forced  convection  but  fell  about  10  percent  above  the  analytical 
curve .  Quite  large  circumferential  variations  in  wall  temperature  were 
observed  in  downflow  as  compared  vfith  those  encountered  in  upflow,  and 
the  fully  developed  Nusselt  numbers  reported  are  based  on  average  wall 
temperatures  determined  by  averaging  the  readings  of  two  diametrically 
opposite  wall  thermocouples  at  each  axial  position.  With  larger  heating 
rates  in  downflow  the  wall  temperature  distributions  strongly  suggested 
a  cell  flow  near  the  bottom.  At  still  larger  heating  rates  the  wall  tem- 
peratur-es  varied  in  a  periodic  way. 


4Tie  information  presented  herein  is  an  abbreviated  form  of  a  thesis 
entitled  "Combined  Forced  and  Free  Convection  in  a  Vertical  Tube"  sub¬ 
mitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy  at  Purdue  University,  1958. 
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imorucTioN 

Description  of  the  Problem 

A  description  of  the  problem  considered  in  this  report  is  perhaps 
best  introduced  by  considering  the  simple  physical  system  illustrated  in 
figure  l(a).  This  sketch  illustrates  a  round  tube  of  finite  length  im¬ 
mersed  in  a  large  body  of  fluid  which  is  at  some  uniform  temperature. 

The  tube  is  positioned  with  its  axis  vertical,  and  the  outer  surface  of 
the  tube  is  perfectly  insulated. 

If  there  is  no  heat  transfer  to  or  work  done  on  the  system,  the 
fluid  remains  at  rest,  the  buoyant  forces  on  each  element  of  fluid  being 
exactly  balanced  by  the  hydrostatic  pressure.  However,  if  the  tube  is 
heated  in  some  way,  the  fluid  within  the  tube  changes  temperature  and 
its  density  decreases. ^  This  reduced  density  causes  an  increase  in  the 
buoyant  force,  and  the  fluid  tends  to  rise  and  pass  out  the  top  of  the 
tube,  where  it  mixes  with  the  large  mass  of  unheated  surrounding  fluid. 
This  action  draws  cold  fluid  into  the  bottom  of  the  tube,  and  hence  a 
pumping  action  occurs .  This  type  of  flow  and  heat  transfer  is  commonly 
called  free  convection. 

If  the  surrounding  fluid  is  replaced  by  the  cooler  piping  system  and 
pump  shown  in  figure  1(b)  and  the  cooling  and  pumping  rates  are  con¬ 
trolled  to  maintain  the  temperatures  ta  and  tg  \inchanged,  conditions 
are  unchanged  inside  the  pipe.  This  situation  is  commonly  called  com¬ 
bined  forced  and  free  convection. 

There  is  one  important  distinction  between  the  two  situations  just 
described.  In  the  first  situation,  for  a  given  heating  rate  in  the  tube, 
there  is  a  unique  flow  rate  established,  which  is  such  that  the  pressure 
drop  through  the  tube  is  exactly  balanced  by  the  net  buoyant  force.  In 
the  second  situation,  on  the  other  hand,  there  is  freedom  to  adjust  the 
pump  flow  rate  to  any  desired  value.  Of  course,  when  the  flow  rate  is 
changed,  the  heat-transfer  rate  must  be  changed  in  order  to  keep  ta 
and  tg  unchanged.  Tiie  buoyant  forces  still  exist,  regardless  of  the 
flow  rate,  but  only  one  flow  rate  corresponds  to  the  free-convection 
situation.  There  is  no  basic  difference  in  the  case  of  "pure  free  con¬ 
vection";  it  is  merely  a  special  case  of  the  so-called  combined  forced 
and  free  convection. 

When  the  flow  rate  is  gradually  increased  in  the  apparatus  shown  in 
figure  l(b),  the  buoyant  force  becomes  less  and  less  important  in  com¬ 
parison  with  the  pumping  force.  The  term  "forced  convection"  is  used  to 
describe  the  limiting  situation  when  the  buoyant  forces  are  negligible. 
Since  all  fluids  change  in  density  when  heated,  pure  forced  convection 


p 

Most  fluids  behave  in  this  manner,  although  there  are  exceptions. 
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is  never  exactly  realized.  One  of  the  goals  of  the  present  study  was  to 
establish  when  the  buoyant,  or  free-convection  effect,  is  negligible. 

There  are  a  great  many  possible  variations  to  the  problem  consid¬ 
ered.  The  tube  walls  may  be  heated  in  different  ways  (two  common  ways 
are  uniform  wall  temperature  and  uniform  wall  heat  flux).  In  addition, 
internal  heat  sources,  due  to  chemical  or  nuclear  reactions,  either  uni¬ 
form  or  nonuniform,  may  be  present  within  the  fluid  inside  the  tube.  The 
tube  diameter  and  length  may  be  varied,  as  well  as  the  kind  of  fluid. 
Finally,  the  cross-sectional  shape  may  be  round,  square,  or  some  other 
shape . 

When  the  cool  fluid  enters  the  bottom  of  the  tube,  it  is  heated  and 
rises.  If  the  tube  is  long  enough,  and  if  the  boundary  conditions  are 
suitable,  the  radial  velocity  and  temperature  profiles  beyond  a  certain 
inlet  length  are  similar  to  the  corresponding  profiles  at  other  cross 
sections  farther  up  the  tube.  This  region  is  known  as  the  region  of 
"fully  developed"  flow  and  heat  transfer.  The  length  of  the  tube  re¬ 
quired  to  establish  similar  profiles  is  called  the  "entrance  length." 

The  problem  considered  herein  is  one  of  combined  forced  and  free 
laminar  convection  in  a  vertical  round  tube.  The  boundary  conditions 
considered  are  uniform  wall  heat  flux.  No  heating  was  present  within 
the  fluid. 

The  transfer  of  heat  by  combined  forced  and  free  convection  in  ver¬ 
tical  channels  (more  generally,  in  channels  in  which  the  flow  is  parallel 
to  the  direction  of  the  body  force,  i.e.,  gravity  or  centrifugal)  is  very 
common,  but  only  in  recent  years  has  much  attention  been  placed  on  deter¬ 
mining  the  characteristics  of  such  systems.  Possible  uses  in  nuclear 
reactor  or  turbine  blade  cooling  are  mainly  responsible  for  recent  inter¬ 
est  in  this  mode  of  heat  transfer. 


Purpose  of  Experiment 

The  purpose  of  the  experiment  was  to  check  the  analysis  of  refer¬ 
ence  1  and  also  to  provide  information  about  the  region  of  applicability, 
that  is,  the  transition  Reynolds  and  Rayleigh  numbers  and  the  entrance 
lengths  required  for  the  flow  to  become  fully  developed.  The  limiting 
case  of  pure  laminar  forced  convection,  including  thermal  entrance  ef¬ 
fects,  was  also  to  be  studied. 

The  specific  case  chosen  for  the  experiment  was  combined  forced  and 
free  laminar  convection  in  a  vertical  round  tube  with  uniform  wall  heat 
flux  and  no  internal  heat  generation.  Provision  was  made  for  both  up- 
flow  with  heating  and  downflov?  with  heating,  that  is,  both  positive  and 
negative  Rayleigh  numbers.  Tne  apparatus  was  designed  to  be  suitable 
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as  well  for  obtaining  data  on  pure  forced  convection,  both  laminar  and 
turbulent.  Previous  data  (refs.  2  and  3)  for  the  laminar  region  were 
subject  to  large  experimental  error. 


SURVEY  OF  THE  LITERATURE 
Theoretical  Investigations 

The  analysis  of  the  case  of  combined  forced  and  free  convection  in 
a  vertical  round  tube  with  iiniform  wall  heat  flux  is  considered  in  ref¬ 
erences  1  and  4  to  8.  None  of  these  analyses  consider  entrance  effects 
or  turbulence  or  transition  to  turbulence .  The  limiting  case  of  pure 
forced  convection  in  thermal  entrance  regions  is  analyzed  in  reference  9. 


Experimental  Investigations 

Some  data  on  combined  forced  and  free  convection  in  vertical  pipes 
with  uniform  wall  heat  flux  are  given  in  references  2,  3,  10,  and  11. 
Reference  11  includes  some  information  on  transition  to  turbulence. 
Since  wall  temperatures  were  not  measured  it  was  not  possible  to  calcu¬ 
late  Nusselt  numbers. 


Comparison  of  Previous  Experiments  with  Analysis 

The  experimental  data  from  other  experiments  are  compared  with  the 
analysis  of  reference  1  in  figure  2.  The  solid  line  represents  the  anal¬ 
ysis.  The  Nusselt  numbers  in  the  range  100  <  Ra  ^  10,000  are  accurately 
approximated  by  the  equation 


Nu  =>  1.40  RaO-28  (l) 

This  relation  is  shown  in  figure  2  by  the  dashed  line.  (Symbols  are  de¬ 
fined  in  appendix  A.) 

The  experimental  points  shown  as  squares  in  figure  2  are  the  data 
of  reference  10  for  high-pressure  water  (2000  and  1500  Ib/sq  in.  abs) 
flowing  upward  in  a  vertical  electrically  heated  tube.  Kie  tube  had  an 
inside  diameter  of  0.1805  inch  and  a  length  of  9  inches  and  thus  an  l/D 
of  50. 

These  data  fall  below  the  analysis.  In  reference  12  these  same  data 
are  compajred  on  a  different  basis  with  other  experimental  data  and  eure 
found  to  fall  low  there  as  well. 
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Two  other  experiments  have  provided  data  for  lower  Rayleigh  num¬ 
bers.  These  experiments  were  designed  for  pure  forced  convective  flow, 
but  the  Rayleigh  numbers  were  high  enough  in  some  cases  to  show  signifi¬ 
cant  free  convective  effects.  The  first  experiment  is  that  of  Gross 
(ref.  2).  These  data  are  shown  as  circles  in  figure  2  and  are  for  water 
in  upflow  in  an  electrically  heated  tube.  Three  tubes  of  0.2244-, 
0.305-,  and  0.432-inch  inside  diameter  were  used. 

Heated  tube  lengths  varied  from  2  to  5  feet.  Local  heat-transfer 
data  were  reported  but  only  data  which  satisfy  the  condition 

>  0.085 

are  plotted  here.^  This  is  the  condition  found  necessary  in  reference 
9  for  the  Nusselt  nximbers  for  pure  forced  convection  larndnar  heat  trans¬ 
fer  to  be  within  5  percent  of  the  fully  developed  value.  This  condition 
is  for  a  fully  developed  entering  velocity  profile  and  a  uniform  wall 
heat  flux  corresponding  to  this  experiment.  Also,  all  data  for  Reynolds 
numbers  greater  than  2000  were  not  plotted.  It  is  seen  that  the  scatter 
in  the  data  of  reference  2  is  quite  large.  The  higher  values  of  Nusselt 
numbers  at  low  and  high  Rayleigh  numbers  are  believed  to  be  caused  by 
bus  bar  conduction  lowering  the  tube  wall  temperatures  at  entrance  and 
exit.  Despite  the  large  scatter,  the  trend  in  the  data  away  from  the 
pure  forced  convection  Nusselt  number  of  4.36  is  largely  accounted  for 
by  the  analytical  curve  shown. 

Two  data  points  at  low  Rayleigh  nujiiber  from  an  experiment  reported 
in  reference  3  are  shown  as  triangles  in  figure  2 .  These  data  are  for 
v/-ater  flowing  up\-mrd  in  an  electrically  heated  vertical  tube.  The  tube 
had  an  inside  diameter  of  0.2244  inch  and  a  heated  length  of  12  feet. 

The  entering  velocity  profile  was  probably  nearly  uniform,  so  the  only 
data  shown  are  those  for  which 


This  is  the  requirement  for  the  Nusselt  numbers  to  be  within  5  percent 
of  the  fully  developed  value  if  the  entering  velocity  profile  is  uniform 
and  the  wall  heat  flux  is  uniform  (ref.  13). 

^For  this  reason  fewer  points  are  shown  in  figure  2  than  are  shown 
in  figure  7  of  reference  1. 


DESCRIPTION  OP  APPARATUS 


Flow  Loop 

A  schematic  diagram  of  the  heat-transfer  loop  is  shown  in  figure  3. 
Distilled  and  deionized  water  was  circulated  as  the  heat-transfer  medium. 
A  constant-head  tank  was  used  to  ensure  a  steady  flow  through  the  test 
section.  This  was  found  to  be  necessary  from  preliminary  tests.  The 
water  was  de-aerated  by  being  heated  to  the  boiling  point  in  the  head 
tank,  and  it  was  then  cooled  to  the  desired  temperature  by  water  cooling 
coils  directly  beneath  the  head  tank.  The  flow  was  down  a  44-foot  in¬ 
sulated  standpipe  to  an  insulation  tank  which  contained  the  heated  test 
section.  The  test  section  ’.^as  surrounded  by  seven  concentric  aluminum- 
foil  radiation  guards.  The  tank  air  pressure,  was  maintained  at  about  1 
micron  of  mercury.  Such  an  arrangement  gave  a  rapid  time  response  and 
nearly  negligible  heat  losses . 

In  upflow  the  fluid  passed  first  through  a  bottom  mixing  box,  where 
the  inlet  bulk  temperature  was  measured,  and  then  up  through  a  3-foot 
insulated  hydrodynamic  approach  section  of  5/l6-inch  bore.  This  section 
had  a  length -to -diameter  ratio  of  115  to  allow  the  velocity  profile  to 
become  fully  developed.  The  water  then  flowed  without  interruption 
through  a  Teflon  insulating  disk  and  into  the  heated  test  section.  The 
Inconel  test  section  also  had  a  5/l6-inch  bore  and  was  3  feet  long 
(l/D  =  115).  It  was  heated  by  passing  low-voltage  alternating  current 
axially  through  the  tube  wall.  Thermocouples  attached  to  the  outer  wall 
allowed  a  determination  of  local  wall  temperature. 

The  heated  water  passed  out  through  a  second  Teflon  insulating  disk 
and  then  to  an  outlet  mixing  box,  where  the  outlet  bulk  temperature  was 
measured.  The  water  left  the  insulation  tank  at  the  top  and  was  cooled 
to  room  temperature  in  a  heat  exchanger.  It  was  then  filtered  and 
passed  through  a  flow  control  needle  valve. 

Tlae  discharge  from  the  flow  control  needle  valve  went  into  two 
glass  burettes,  where  the  flow  rate  was  determined.  This  measurement 
was  accomplished  by  determining  the  time  required  to  fill  the  burettes 
when  the  outlet  valve  was  closed.  The  time  was  measured  automatically 
by  a  combination  of  electronic  relays  and  an  electronic  timer.  Figure 
3  illustrates  this  system.  A  rotameter  was  also  used  to  give  a  rough 
indication  of  the  flow  rate. 

The  water  leaving  the  burette  went  to  a  sump  tank  eind  then  was 
pumped  back  up  to  an  overflow  tank  by  a  sump  pump.  A  second  circulating 
pump  transferred  the  water  from  the  overflow  tank  through  a  heated  pipe 
to  the  head  tank.  The  purpose  of  the  heater  was  to  heat  the  water  to 
near  the  boiling  point  and  thus  de-aerate  it.  The  overflow  from  the 
head  tank  returned  to  the  overflow  tank  through  a  cooler.  The 
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temperature  of  the  water  in  the  head  tank  was  held  near  the  boiling 
point  by  a  thermoswitch-heater  combination. 

A  bypass  mixed-bed  ion  exchanger  not  shown  in  figure  3  was  con¬ 
nected  between  the  bottom  of  the  standpipe  and  the  sump  tank.  A  con¬ 
stant  flow  through  the  ion  exchanger  maintained  the  high  purity  of  the 
distilled  water. 

Details  of  the  various  components  of  the  flow  loop  are  given  in 
reference  3. 

For  downflow  the  water  entered  the  insulation  tank  at  the  top  and 
left  at  the  bottom.  The  rest  of  the  flow  loop  was  unchanged.  The 
starting  length  for  downflow  was  considerably  shorter,  being  about  13 
diameters . 

Several  modifications  were  made  to  the  apparatus  after  17  runs  had 
been  made.  One  change  was  the  addition  of  a  temperature  regulating  bath 
to  control  water  inlet  temperature.  This  bath  reduced  inlet  temperature 
fluctuations  from  1°  F  to  less  than  0.2°  F.  Another  modification  was 
the  installation  of  an  immersion  heater  in  the  head  tank  to  replace  a 
wrap-around  heater  on  the  inlet  line.  This  change  allowed  more  reliable 
operation  and  a  steadier  head  tank  temperature .  A  third  modification 
was  an  improvement  in  the  vibration  isolation  of  the  insulation  tank. 
Vibrations  from  the  vacuuju  pump  had  been  felt  in  the  insulation  tank  and 
there  was  a  possibility  that  this  could  cause  an  early  transition  to 
turbulent  flow. 


Test  Section  and  Approach  Section 

The  test  section  was  a  resistance-heated  Inconel  tube  with  a 
0.3130-inch  bore  and  a  0.0312-inch  wall  thickness.  The  length  was 
36.030  inches,  and  the  length-to -diameter  ratio  was  therefore  about  115. 
A  sketch  of  the  test  section  and  approach  section  is  shown  in  figure  4, 
and  a  photograph  of  the  assembly  is  given  in  figure  5.  The  approach 
section  in  downflow  had  a  length-to-diameter  ratio  of  13. 

Fifty-eight  iron-constantan  thermocouples  were  spot-welded  to  the 
outside  of  the  test  section  in  29  diametrically  opposed  pairs.  A  detail 
of  the  method  of  attachment  is  shown  in  figure  4.  A  listing  of  their 
exact  locations  is  given  in  reference  5.  The  axial  spacing  of  the  ther¬ 
mocouples  near  each  end  was  small  in  order  to  allow  an  accurate  determi¬ 
nation  of  entrance  region  wall  temperatures  for  both  upflow  and  downflow 
Axially  adjacent  pairs  of  thermocouples  were  positioned  90°  apart.  This 
was  done  to  facilitate  detection  of  any  asymmetry  in  wall  temperature 
distribution. 
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Circular  copper  flanges  were  silver  soldered  to  each  end  of  the  test 
section.  The  power  leads  were  connected  to  the  outer  circumference  of 
these  flanges.  Potential  leads  were  soldered  to  each  flange.  These 
flanges  also  contained  thermocouples  which  were  used  to  determine  heat 
losses  as  is  described  in  reference  5.  The  test  section  was  connected 
to  and  insulated  from  the  approach  and  exit  sections  by  Teflon  disk 
spacers . 

A  detail  drawing  of  the  approach  section  is  given  in  figure  4.  The 
approach  section  was  identical  to  the  test  section  except  that  it  was 
not  heated  and  had  fewer  thermocouples  attached  to  it.  The  approach  sec¬ 
tion,  Teflon  spacers,  test  section,  and  exit  section  were  all  very  care¬ 
fully  alined,  doweled,  and  bolted  in  order  to  ensure  that  the  flow  pat¬ 
tern  would  not  be  disturbed  when  the  water  passed  from  one  section  to 
another . 


Test-Section  Insulation 

The  small  heat-transfer  coefficients  encountered  in  laminar  flow  re¬ 
quire  an  insulation  superior  to  glass  or  quartz  wool.  Two  possibilities 
which  seemed  promising  were  a.  vacuujn  bottle  or  a  powdered  insulation 
under  vacuum.  The  first  approach  was  followed,  because  it  was  felt  that 
the  large  time  constants  involved  with  a  low-thermal-diffusivity  material 
would  make  operation  difficult. 

Seven  concentric  cylinders  of  aluminum  foil  were  mounted  between  the 
test  section  and  the  insulation  tank  to  act  as  radiation  shields . 

The  radial  heat  loss  was  measured  with  the  test  section  drained  and 
plugged.  A  plot  of  heat  loss  per  foot  is  shown  in  figure  6.  The  loss  is 
nearly  a  linear  function  of  temperature  difference.  It  is  possible  to 
compute  a  "thermal  conductivity"  of  the  insulation  tank.  It  is  about 
0.0046  Btu/(hr)(ft)(°F)  compared  with  0.023  Btu/(hr)(ft)(OF)  for  glass 
wool.  Thus  it  is  about  five  times  as  effective.  The  measurements  were 
all  made  at  the  tank  operating  pressure  of  0.6  to  1  micron  of  mercury. 
Tliis  radial  heat  loss  is  nearly  negligible  for  most  nms  but  was  included 
in  all  calculations. 


Power  Supply 

The  test  section  was  resistance  heated  by  the  power  supply  shown 
schematically  in  figure  7 .  Power  was  supplied  from  a  230-volt  supply, 
and  the  current  passed  through  a  5 -kilovolt -ampere  electronic  voltage 
regulator  with  a  line  voltage  regulation  accuracy  of  ±0.1  percent.  The 
output  then  passed  to  a  variable  autotrans former  and  then  through  two  2- 
kilovo It -ampere  stepdown  transformers  in  parallel,  which  reduced  the 


voltage  to  a  maximum  of  11  volts.  From  the  transformers  the  current 
flowed  directly  to  the  power  lead  connections  on  top  of  the  Insulation 
tank. 
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Power  Measuring  Instruments 

Power  was  measured  by  a  multiplied-deflection  a-c  potentiometer  sim¬ 
ilar  to  the  one  described  in  reference  14.  This  method  has  the  advantage 
that  the  accuracy  is  limited  only  by  the  full-scale  accuracy  of  the  me¬ 
ters  and  is  independent  of  power  level.  This  was  highly  desirable  be¬ 
cause  power  ranged  from  0.17  watt  on  the  heat  loss  measurements  to  4000 
watts  on  the  higher  power  heat-transfer  runs.  An  a-c  potentiometer  was 
of  value  also,  because  the  low  voltage  drops  across  the  test  section 
made  it  difficult  to  use  a  wattmeter  directly.  The  voltmeter  and  ammeter 
were  calibrated  in  an  NASA  instruments  laboratory  and  correction  charts 
were  used  with  them.  The  wattmeter  was  calibrated  with  a  resistor  whose 
d-c  resistance  was  accurately  measured  on  a  Kelvin  bridge.  The  wattmeter 
readings  were  compared  with  of  the  test  section,  and  they  always 

agreed  to  within  2  percent.  The  power  factor  was  very  nearly  unity  at 
all  times . 


Temperature  Measuring  System 

Iron-constantan  thermocouples  were  used  throughout  for  temperature 
measurement.  Test-section  and  mixing-box  thermocouples  were  carefully 
calibrated.  Details  of  t;he  calibration  procedure  are  given  in  refer¬ 
ence  5.  Tlie  tables  of  ITBS  circular  561  were  used  for  uncalibrated  ther¬ 
mocouples  . 

A  circuit  diagram  i'or  the  thermocouples  is  shown  in  figure  8.  Test- 
section  wall  thermocouples  and  the  outlet  mixing-box  thermocouple  were 
always  read  with  respect  to  the  inlet  mixing-box  thermocouple.  Since 
both  upflow  and  downflow  runs  were  made,  provision  was  made  for  a  re¬ 
versing  switch  to  interchange  roles  of  the  mixing  boxes .  Each  mixing- 
box  junction  could  be  read  with  respect  to  an  ice  junction  in  order  to 
establish  the  absolute  temperature  level  of  the  system.  A  self-balancing 
potentiometer  was  used  to  obtain  rough  readings  when  the  system  was  ap¬ 
proaching  steady  state. 

A  Rubicon  potentiometer.  Model  B,  was  used  to  measure  thermocouple 
voltages .  An  electronic  null  indicator,  with  a  voltage  sensitivity  of 
1  microvolt  per  small  division,  served  in  place  of  a  conventional  light 
beam  galvanometer  to  indicate  the  matching  of  potentiometer  and  thermo¬ 
couple  volteiges . 
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All  thermocouple  wires  passed  to  a  junction  box.  The  thermocouple 
vires  on  the  test  section  and  in  the  mixing  boxes  were  36  gage  and  were 
all  cut  from  the  same  spools  of  wire .  At  the  junction  box  the  wires  were 
increased  in  size  to  make  the  leads  more  durable.  Iron-constantan  wire 
was  used  throughout  the  thermocouple  circuits.  From  the  junction  box  the 
leads  vent  to  a  switch  box,  which  housed  all  switches .  Leads  from  the 
switch  box  passed  to  the  self-balancing  potentiometer  and  the  Rubicon 
potentiometer . 

The  junction  box  and  switch  box  were  of  identical  construction.  All 
thermocouple  junctions  were  contained  in  a  copper  box.  This  copper  box 
was  surrounded  by  several  inches  of  foamed-plastic  insulation,  which  in 
turn  was  contained  in  a  wooden  box.  The  purpose  of  this  type  of  con¬ 
struction  was  to  produce  a  uniform  temperature  at  the  junction  of  each 
pair  of  thermocouple  wires .  Even  though  the  box  temperature  might  chajige 
because  of  changes  in  ambient  conditions,  both  thermocouple  junctions 
would  change  in  the  same  way  and  thus  prevent  euiy  stray  voltages  from 
being  developed. 


EXPERIMENTAL  PROCEDURE 

The  experimental  data  were  taken  in  the  following  manner.  Several 
runs  at  different  heating  rates  were  usually  taken  at  one  flow  rate. 

When  a  new  flow  rate  was  set  by  adjusting  the  flow  control  needle  valve, 
it  was  generally  necessary  to  allow  the  apparatus  to  run  overnight  in 
order  that  steady  conditions  would  be  established.  When  the  flow  was 
steady  and  the  inlet  water  temperature  variations  were  tolerable,  the 
power  was  turned  on  and  set  at  the  desired  value.  Usually  a  half  hour 
was  sufficient  to  attain  equilibrium.  This  time  was  adequate  because  of 
the  fast  time  response  of  the  thin  aluminum  radiation  guards . 

Readings  were  taken  of  volts,  amperes,  watts,  voltage  divider  ratio, 
current  transformer  ratio,  water  inlet  temperature  above  ice  point,  water 
outlet  temperature  above  ice  point,  eind  bulk  temperature  rise  before, 
half  way  through,  and  after  the  reading  of  wall  temperature  thermocouples 
As  many  readings  of  the  time  to  fill  burettes  were  taken  as  was  possible 
during  the  run.  Besides  wall  temperature  readings,  data  were  taken  on 
room  temperature,  flange  temperatures,  and  approach  section  temperatures. 

At  the  completion  of  a  run  the  power  level  was  readjusted  to  obtain 
a  different  Rayleigh  number. 


EXPERIMENTAL  RESULTS 

The  important  results  of  the  experiment  are  presented  and  discussed 
in  this  section,  and  a  discussion  of  accuracy  and  anomalies  of  the 
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experiment  is  given  in  appendix  B.  The  experimental  results  are  tabu¬ 
lated  in  reference  5. 


Results  for  Heating  with  Upflow 

Fully  developed  heat  transfer.  -  The  primary  objective  of  this  ex¬ 
periment  was  to  verify  the  analysis  which  was  made  for  the  condition  of 
fully  developed  heat  transfer  and  constant  fluid  properties .  The  data 
which  are  believed  to  be  fully  developed  and  which  were  steady  (nonfluc¬ 
tuating)  are  shown  in  figure  9.  The  solid  line  is  the  predicted  curve 
of  the  analysis,  which  is  shown  in  figure  2  also.  It  is  apparent  that 
these  data  confirm  the  analysis  over  the  range  presented.  The  data  at 
high  values  of  Rayleigh  number  show  some  scatter  and  also  seem  to  be 
slightly  below  the  analytical  curve.  Physical  properties  were  evaluated 
at  the  film  temperature;  some  of  the  scatter  may  be  due  to  physical  prop¬ 
erty  variations  with  temperature.  However,  it  is  believed  that  these 
data  are  low  at  high  Rayleigh  numbers  because  of  a  peculiar  entrance  ef¬ 
fect.  This  is  discussed  later  when  the  thermal  entrance  region  results 
are  presented. 

Another  observation  to  be  made  from  figure  9  is  that  it  is  very  dif¬ 
ficult  to  confirm  the  pure  forced  convection  fully  developed  Nusselt  num¬ 
ber  of  48/ll  (4.36)  experimentally  because  of  the  influence  of  free  con¬ 
vection  effects .  The  lowest  Rayleigh  number  attainable  with  reasonable 
accuracy  was  about  25.  This  gives  a  Nusselt  number  of  4.62,  about  6  per¬ 
cent  above  the  zero  Rayleigh  number  value. 

Tliermal  entr^ce  region  data.  -  Having  confirmed  the  fully  developed 
analysis,  a  question  arises  as  to  what  length  is  required  to  produce  the 
fully  developed  condition.  An  investigation  of  this  thermal  entrance 
length  may  be  made  by  means  of  the  experiment. 

The  pure  forced  convection  (zero  Rayleigh  number)  analysis  of  ref¬ 
erence  9  provides  a  basis  for  comparison  of  thermal  entrance  region  data. 
The  Nusselt  numbers  very  near  the  heated  entrance  will  be  independent  of 
the  heating  rate,  because  the  Nusselt  number  is  determined  only  by  the 
flow  pattern.  The  velocity  profile  will  not  be  distorted  from  its  en¬ 
tering  parabolic  shape  until  a  temperature  pattern  is  created  in  the 
fluid.  This  temperature  pattern  then  creates  radial  density  differences, 
which,  in  turn,  distort  the  parabolic  velocity  profile  and  thus  change 
the  Nusselt  munbers  from  the  pure  forced  convection  values.  A  certain 
heated  length  is  required  for  this  profile  to  be  distorted. 

The  pure  forced  convection  analysis  results  in  a  single  relation 
between  Nusselt  number  and  the  parameter  is  reasonable,  then, 

to  present  the  data  in  terras  of  these  variables,  with  Rayleigh  niamber  as 
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a  parameter.  Such  a  comparison  is  given  in  figure  10.  The  solid  curve 
is  the  zero  Rayleigh  number  analysis  (ref.  9),  and  the  data  shown  are  for 
a  very  low  Rayleigh  number  run  and  for  a  typical  high  Rayleigh  number  run. 
Numbers  beside  the  data  indicate  local  Rayleigh  numbers . 

The  very  low  Rayleigh  number  run  agrees  well  with  the  pure  forced 
convection  analysis  of  reference  9.  These  data  are  typical  of  all  low 
Rayleigh  number  rvins. 

The  high  Rayleigh  number  run  agrees  with  the  zero  Rayleigh  number 
analysis  very  near  the  entrance,  but  soon  breaks  away,  goes  through  a 
minimum,  and  rises  again.  The  dashed  line  is  the  fully  developed  com¬ 
bined  forced  and  free  convection  analysis  (ref.  l)  applied  locally.  The 
Rayleigh  number  increases  toward  the  pipe  exit  because  of  physical  prop¬ 
erty  variations.  The  agreement  between  the  data  and  the  dashed  line  in¬ 
dicates  that  the  concept  of  a  "locally  fully  developed"  heat  transfer  is 
a  valid  one.  Again,  this  run  is  typical  of  many  such  runs  taken  and  tab¬ 
ulated  in  tables  II  and  III  of  reference  5 . 

The  highest  Rayleigh  number  run  in  which  entrance  data  were  obtained 
is  shown  on  the  same  type  of  plot  in  figure  11.  Again  the  analysis 
(ref.  l)  applied  locally  is  shown  as  a  dashed  line.  In  this  run  the  data 
fall  below  the  analysis  and  then  gradually  climb  to  approach  the  analyti¬ 
cal  curve.  The  effect  is  most  pronounced  in  this  highest  Rayleigh  nxim- 
ber  run.  This  effect  was  observed  to  a  lesser  extent,  in  lower  Rayleigh 
number  rxins  as  well.  It  may  be  due  to  experimental  error,  but  this  is 
considered  unlikely.  At  high  values  of  Rayleigh  number  an  "undershoot" 
occurs,  because  the  large  velocity  profile  distortion  involved  requires 
a  considerable  length  in  which  to  develop. 

This  argument  and  discussion  are  based  on  somewhat  meager  data,  but 
they  are  believed  to  explain  the  low  Nusselt  numbers  observed  at  high 
Rayleigh  numbers  in  figure  9.  These  particular  data  were  taken  in  the 
vicinity  of  the  transition  points,  discussed  later,  and  complete  entrance 
data  were  not  obtained.  It  is  likely  that  some  of  these  data  are  in  the 
"undershoot"  region  and  hence  would  fall  below  the  analysis. 

For  pure  forced  convection,  it  is  shown  in  reference  9  that  the 
value  of  (2x/D)/RePr  at  which  the  Nusselt  number  came  to  within  5  percent 
of  its  fully  developed  value  was  0.085. 

The  entrance  length  for  the  high  Rayleigh  number  run  (run  (U14)) 
shown  in  figure  10  was  taken  to  be 
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and  for  run  (U24)(fig.  ll)  this  value  was  taken  to  he  about  0.075  where 
the  data  become  asymptotic  to  the  analysis,  applied  locally.  Values  of 
(2x/D)/RePr  at  which  the  data  became  parallel  to  the  analytical  curve 
are  plotted  in  figure  12.  The  scatter  in  the  data  is  very  large  because 
of  the  difficulty  in  accurately  estimating  when  the  fully  developed  con¬ 
dition  is  reached.  However,  the  trend  definitely  indicates  a  decrease 
in  thermal  entrance  length  with  increase  in  Rayleigh  number  followed  by 
an  increase  at  large  Rayleigh  numbers. 

Transition  data.  -  A  transition  from  steady  laminar  flow  to  a  slow, 
apparently  random,  eddy  flow  was  observed  on  some  runs.  This  was  evi¬ 
denced  by  wall  temperature  fluctuations  appearing  on  the  upper  portions 
of  the  heated  tube.  Temperature  records  taken  with  a  photoelectric  re¬ 
corder  are  shown  in  figure  13  for  a  typical  run  in  which  transition  oc¬ 
curred.  With  an  increase  in  heating  rate,  at  a  constant  flow  rate,  the 
transition  point  would  travel  down  the  tube  toward  the  entrance.  If  the 
flow  rate  were  increased,  keeping  the  Rayleigh  number  constant,  the  tran¬ 
sition  point  would  travel  up  and  eventually  out  the  tube  exit. 

Special  runs  were  made  in  order  to  gather  data  on  this  transition. 

A  tentative  correlation  of  these  transition  data  is  shown  in  figure  14. 

KbPi* 

Plotted  here  against  ^x/d  ^  special  Rayleigh  number  defined  by 

P^Pg(Vw 
Han  =  - 


Pr 


(2) 


Physical  properties  are  evaluated  at  the  film  temperature.  The  dashed 
line  is  given  by  the  equation 


Since  thermocouples  were  located  only  at  2 -inch  intervals  in  most  cases, 
it  was  difficult  to  locate  the  transition  x/D  exactly.  In  addition, 
the  transition  point  was  difficult  to  locate  precisely,  because  the  wall 
temperature  fluctuations  were  damped  by  the  heat  capacity  of  the  tube 
wall.  The  correlation  is  probably  adequate  for  most  applications,  how¬ 
ever  . 


Since  transition  is  likely  to  be  affected  by  such  things  as  vibra¬ 
tion  level,  tube  nonuniformities,  the  tube  being  not  exactly  "vertical, 
and  so  forth,  it  is  recommended  that  additional  data  be  obtained  in  an¬ 
other  apparatus . 

Transition  data  recently  obtained  are  reported  in  reference  11.  The 
transition  point  was  detected  by  injecting  dye  into  the  water  upstream 
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and  observing  when  the  dye  stream  began  to  fluctuate  at  the  exit  of  a 
transparent  tube.  These  data  are  shown  in  figure  14,  Since  observations 
were  made  at  the  exit  only,  no  x/D  variation  could  be  made.  The  tran¬ 
sition  point  for  stream  fluctuation  occurs  at  lower  values  of  Rap  than 
for  wall  temperature  fluctuations.  This  would  be  expected,  since  the 
stream  fluctuations  would  have  to  build  up  to  some  magnitude  in  order  to 
transport  heat  radially  from  the  wall. 

Selection  of  data.  -  A  discussion  of  the  selection  of  upflow  data 
is  perhaps  best  introduced  by  referring  to  two  representative  plots  of 
temperature  against  x/D. 

Figure  15(a)  shows  a  representative  wall  temperature  plot  for  a  low 
Pec let  number  (RePr)  run.  Generally  there  were  two  thermocouples  at  each 
x/D,  located  diametrically  opposite.  In  some  cases,  however,  only  one 
thermocouple  was  operative.  The  dashed  line  shown  is  a  curve  drawn 
through  the  points  obtained  by  averaging  temperatures  only  at  axial  loca¬ 
tions  where  both  thermocouples  were  operative.  The  circles  represent  in¬ 
dividual  thermocouple  temperatures.  Several  observations  can  be  made 
from  this  figure.  First,  a  very  smooth  curve  can  be  drawn  through  the 
averaged  data.  Second,  an  asymmetric  wall  temperature  distribution  is 
apparent  near  the  tube  exit.  Third,  a  slight  rise  in  temperature  exists 
in  the  first  several  diameters. 

The  fact  that  a  smooth  curve  can  be  drawn  through  the  averaged  data 
was  used  to  reject  data  points  obtained  when  only  one  thermocouple  ex¬ 
isted  at  a  location,  since  these  points  generally  did  not  fall  on  the 
smooth  curve  because  of  the  asymmetry  in  wall  temperature.  Tlie  asymmetry 
in  wall  temperature  was  more  severe  at  higher  powers.  This  could  be 
caused  by  a  nonuniform  tube  wall  causing  nonuniform  heating,  or  it  could 
be  caused  by  an  asymmetric  flow  pattern  set  up  by  free  convection  effects. 
It  is  believed  to  be  principally  the  latter,  because  in  downflow  the  asym¬ 
metry,  which  had  been  near  the  top  in  upflow,  now  appeared  near  the  bot¬ 
tom.  This  is  discussed  in  the  next  section.  The  theory  predicts  that 
the  wall  temperatures  should  approach  the  fluid  bulk  temperature  at  the 
entrance.  This  does  not  happen  in  the  experiment;  in  fact,  there  is  a 
slight  temperature  rise  very  near  the  entrance.  This  is  caused  by 
joulean  heat  from  the  bottom  bus  bar  being  conducted  into  the  bottom 
flange.  The  first  few  temperatures  were  inaccurate  because  of  this. 

Figure  15(b)  shows  a  representative  wall  temperature  plot  for  a  high 
Peclet  number  rvin.  The  small  rise  in  temperature  near  the  entrance  is 
still  present.  The  asymmetry  in  wall  temperature  is  now  more  severe  at 
the  higher  power.  The  striking  thing  about  this  plot,  however,  is  the 
dip  in  wall  temperature  near  the  exit  of  the  tube.  Tliis  dip  was  present 
in  all  runs  with  Peclet  numbers  of  1500  or  more.  The  dips  had  a  minimum 
at  a  constant  x/D  of  about  95  (about  6  in.  from  exit)  which  did  not 
seem  to  shift  under  any  circumstances.  The  explanation  for  this  dip  is 


15 


unknown.  It  has  been  established  that  it  is  not  (l)  due  to  heat  conduc¬ 
tion  through  a  bundle  of  thermocouple  leads  which  leaves  at  this  pointy 
(2)  due  to  a  heat  leak  in  the  insulation  at  this  point,  (3)  caused  by  a 
nonuniform  tube  in  this  region,  or  (4)  a  false  temperature  caused  by  a-c 
effects  on  the  thermocouple  circuits.  No  dip  occurred  in  any  downflow 
runs.  Data  in  this  region  are  not  presented  in  the  results  (although 
they  are  tabulated  in  ref.  5)  because  the  reason  for  this  behavior  is 
not  known. 


Results  for  Heating  with  Downflow 

The  downflow  runs  were  characterized  by  an  asymmetry  in  wall  temper¬ 
ature  which  became  more  severe  as  the  heating  rate  was  increased.  For 
the  highest  heating  rates,  the  flow  became  unsteady,  and  the  wall  temper¬ 
atures  in  the  lower  portion  of  the  heated  tube  oscillated  in  a  periodic 
way. 


Fully  developed  heat  transfer.  -  The  data  which  appeared  to  be  fully 
developed  are  presented  in  figure  16  as  a  plot  of  Nusselt  number  against 
Rayleigh  number.  The  line  with  which  the  data  are  to  be  compared  is  the 
solid  line  labeled  negative  Rayleigh  number.  This  line  is  derived  in 
reference  5.  Tlie  two  dashed  lines,  for  zero  Rayleigh  number  and  for  pos¬ 
itive  Rayleigh  number,  are  included  to  provide  a  comparison  with  previ¬ 
ously  discussed  upflow  results. 

It  is  seen  that  most  of  the  data  fall  about  10  percent  above  the 
predicted  curve.  Each  point  shown  is  based  on  a  wall  temperature  obtained 
by  averaging  the  temperatures  of  two  thermocouples  placed  diametrically 
opposite.  These  temperature  differences  were  as  high  as  7°  F  in  some 
cases  when  the  axial  temperature  rise  over  the  length  of  the  tube  was 
only  2'tP  F.  It  is  believed  that  this  asymmetry  in  wall  temperature  is 
caused  by  an  asymmetry  in  flow  pattern.  It  is  not  surprising,  therefore, 
that  the  data  do  not  eigree  exactly  with  the  analysis . 

Even  though  most  of  the  data  fall  above  the  analysis,  the  trend  in 
the  data  is  correctly  predicted.  All  the  Nusselt  numbers  are  below  the 
pure  forced  convection  value  of  4.36,  and  considerably  below  the  data 
which  have  exactly  the  same  heating  and  flow  conditions  but  which  have 
the  direction  of  flow  reversed. 

Thermal  entrance  region  data.  -  A  typical  plot  of  Nusselt  number 
against  the  parameter  (2x/D)/RePr  is  shown  in  figure  17.  The  data  do  not 
reach  the  fully  developed  analysis  (dashed  line)  for  reasons  already  dis¬ 
cussed.  The  entrance  lengths  appear  to  be  longer  than  for  pure  forced 
convection,  although  it  is  difficult  to  estimate  entrance  lengths  accu¬ 
rately.  Estimation  is  difficult  raeunly  because  the  local  Rayleigh  number 
decreases  (negatively)  with  length  and  so  the  "local  fully  developed" 
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Nusselt  numbers  decrease  with  length.  The  range  of  Rayleigh  numbers  was 
quite  small.  A  rough  estimate  for  entrance  length  is 


if  the  Rayleigh  number  is  not  below  -120. 

Closer  to  the  entrance,  the  Nusselt  numbers  agree  quite  well  with 
the  pure  forced  convection  analysis  (as  they  did  for  positive  Rayleigh 
numbers).  The  two  low  points,  nearest  the  entrance,  are  believed  to  be 
caused  by  bus  bar  conduction  errors . 

The  approach  section  length-to-diameter  ratio  for  downflow  was  about 
13.  The  length-to-diameter  ratio  required  for  fully  developed  flow  is 
about  0.0575  Re  (ref.  15).  Most  of  the  runs  were  made  with  fully  devel¬ 
oped  flow  at  the  entrance  to  the  heated  section.  The  highest  Reynolds 
number  used  was  330,  which  would  require  an  L/D  of  19  for  fully  devel¬ 
oped  flow. 

Axially  asymmetric  and  periodic  flow.  -  The  wall  temperature  was 
found  to  vary  around  the  circumference  of  the  tube  foi  all  downflow  runs. 
For  a  given  run,  the  variations  became  worse  from  inlet  to  outlet  (top  to 
bottom),  and  for  different  runs,  the  variations  became  more  severe  with 
increasing  heating  rates.  For  low  power  runs,  the  averaging  of  diamet¬ 
rically  opposite  thermocouple  readings  gave  a  smooth  curve  from  top  to 
bottom  (fig.  18(a)),  whereas  for  higher  power  runs  an  apparent  discon¬ 
tinuity  in  this  "average"  wall  temperature  curve  is  indicated  (fig. 
18(b)).  The  variation  across  the  tube  below  this  discontinuity  became 
very  large  in  some  cases,  being  as  high  as  20°  F.  One  might  speculate 
that  above  this  discontinuity  the  flow  is  always  down,  but  is  asymmet¬ 
ric;  below  the  discontinuity  a  cell  flow  may  be  set  up,  with  a  large 
crossflow  at  the  position  of  the  discontinuity.  Since  it  was  not  pos¬ 
sible  to  observe  the  flow,  however,  one  can  only  speculate. 

With  a  still  greater  increase  in  heating  rate,  the  wall  temperatures 
became  unsteady  in  the  lower  portions  of  the  tube.  The  variations  were 
nearly  periodic  in  nature .  The  amplitude  of  the  variations  was  different 
at  different  points  on  the  tube,  but  the  frequency  was  the  same  every¬ 
where  .  Figure  19  gives  a  record  of  the  fluctuations  observed  for  one 
particular  run.^  The  variations  were  more  regular  in  the  upper  part  of 
the  region  than  in  the  lower  part.  The  frequency  was  about  0.8  cycle  per 
minute  for  this  example  and  was  about  0.6  cycle  per  minute  for  one  other 
run  recorded. 


^In  figure  19,  the  indicated  position  of  the  recordings  with  re¬ 
spect  to  circumferential  location  is  not  significant.  The  leaders  are 
used  only  to  indicate  Eocial  position. 
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Selection  of  data.  -  In  downflow  the  circumferential  wall  tempera¬ 
ture  variations  were  much  more  severe  than  they  were  in  upflow.  For  low 
heating  rates  the  averaging  of  data  for  pairs  of  thermocouples  at  each 
axial  position  resulted  in  a  smooth  average  wall  temperature  curve,  as 
indicated  by  figure  18(a).  In  that  figure,  symbols  of  the  same  kind  are 
all  on  the  same  side  of  the  tube  6ind  a  smooth  curve  could  be  drawn 
through  each  set  of  symbols.  The  temperature  differences  across  the  tube 
are  as  high  as  4°  F  for  this  particular  run. 

A  different  wall  temperature  behavior  was  observed  at  higher  heating 
rates,  as  shown  in  figure  18(b).  The  points  can  be  averaged  to  give  a 
smooth  curve  up  to  a  certain  position,  but  then  very  large  circumferen¬ 
tial  wall  temperature  differences  exist  for  the  rest  of  the  tube,  the 
differences  being  as  high  as  9°  F.  These  temperatures  are  steady  in  time. 
Tlie  temperature  pattern  would  indicate  a  cell  flow.  With  a  still  further 
increase  in  power,  the  periodic  fluctuations  discussed  previously  ap¬ 
peared. 

Since  the  flov/  is  now  down  instead  of  up,  the  exit  of  the  tube  be¬ 
comes  the  entrance.  The  asymmetries  present  in  upflow  do  not  exist  at 
the  top  of  the  tube,  nor  was  a  dip  in  wall  temperature  ever  observed 
there . 

Only  wall  temperatures  which  fit  the  smooth  curve  before  the  sudden 
change  in  wall  temperature  distribution  were  used  in  plotting  results  to 
be  compared  vdth  the  analysis.  All  data  are  presented  in  tables  V  and 
VI  of  reference  5.  The  temperature  of  each  individual  wall  thermocouple 
is  reported  in  table  VI  of  reference  5  along  with  the  average  wall  tem¬ 
perature  at  each  location. 


CONCLUDIITG  REMARKS 

The  fully  developed  solutions  of  the  analysis  are  all  in  the  nature 
of  asymptotic  solutions  and  are  expected  to  apply  only  far  from  the  en¬ 
trance.  The  results  of  the  experiment  indicate  that,  for  positive 
Rayleigh  numbers,  these  solutions  do  apply,  and  the  concept  of  locally 
fully  developed  results  is  established.  It  is  also  concluded  that  the 
boundary  conditions  associated  with  these  solutions  are  physically 
realizable . 

The  results  of  the  analysis  were  found  to  apply  only  for  values  of 
negative  Rayleigh  nunber  which  are  small;  for  large  negative  Rayleigh 
numbers  axial  asymmetry  and  unsteady  flow  were  observed,  and  these  con¬ 
ditions  make  the  analysis  inapplicable. 

The  thermal  entrance  region  analysis  of  Siegel,  Sparrow,  and  Hallman 
for  pure  forced  convection  was  experimentally  confirmed.  However,  the 
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fully  developed  Nusselt  niimber  of  48/ll,  predicted  for  this  situation,  is 
difficult  to  verify  experimentally  because  of  the  influence  of  free  con¬ 
vection  effects . 

The  transition  observed  during  heating  with  upflow  would  lead  to  the 
conclusion  that  similar  phenomena  are  likely  to  be  encountered  in  other 
similar  types  of  flows.  This  would  be  expected  to  occur,  for  example,  in 
upflow  of  a  heat-generating  fluid  in  an  insulated  pipe.  It  is  recom¬ 
mended  that  further  experiments  be  run  to  establish  the  region  of  appli¬ 
cability  of  those  portions  of  the  analysis  made  in  the  Hallman  thesis 
which  were  not  investigated  by  the  present  experiment. 


Lewis  Reseajrch  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  August  24,  1961 
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APPENDIX  A 
SYMBOIfi 

A  axial  temperature  gradient  in  fluid,  St/5x,  °F/ft 
Cp  specific  heat  of  fluid  at  constant  pressure,  Btu/(l'b)(OF) 

D  tube  inside  diameter,  ft 

g  acceleration  due  to  gravity,  ft/ sec^ 

h  heat-transfer  coefficient,  q/^m,  Btu/(sec)(sq  ft)(OF) 

I  test-section  current,  amp 

k  thermal  conductivity  of  fluid,  Btu/(sec)(sq  ft)(°F/ft) 

kf  thermal  conductivity  evaluated  at  film  temperature 

L  length  of  heated  section,  ft 

Nu  Nusselt  number,  hD/k,  dimensionless 

Pr  Prandtl  number,  ^lgCp/k,  dimensionless 

q  wall  heat-flux  density,  Btu/(sec)(sq  ft) 

q^  net  electrical  input  to  fluid,  Btu/hr 

qth  thermal  energy  rise  of  fluid  in  passing  through  heated  test  sec¬ 
tion,  Btu/hr 

R  test-section  resistance,  ohms 

Ra  Rayleigh  number,  p2pg2cpD^A/l6nk,  dimensionless 

Raj)  Rayleigh  number  based  on  dieuneter  of  tube,  p^Pg^CpD^Cty  - 
dimensionless 

Re  Reynolds  number,  u^D/v,  dimensionless 

t  static  fluid  temperature,  °F 

tg,  ambient  temperature,  °F 

tg  exit  temperature,  °F 
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"tin 

X 

3 


mixed-mean  or  bulk  fluid  temperature,  °F 

inside  tube  wall  temperature,  °F 

mean  velocity  of  fluid  in  tube,  ft/sec 

distance  measured  along  axis  of  tube  upward 
length,  ft 

thermal  coefficient  of  volume  expansivity,  - 


from  start  of  heated 

ft/°F)/cu  ft 


0,,^  mixed-mean  to  wall  temperature  difference,  °F 
p  dynamic  viscosity  of  fluid,  (lb )( sec  )/sq_  ft 
V  kinematic  viscosity  of  fluid,  sq  ft/sec 
p  mass  density  of  fluid,  ( lb )(sec2 )/ft^ 

Subscripts ; 

f  fiLn  temperature 


fd 


fully  developed 
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APPENDIX  B 
ACCURACY  OP  RESULTS 
Error  Analysis 

A  detailed  error  analysis  was  made  for  one  set  of  thermocouples  in 
one  run  (run  (U30)),  and  the  result  for  Nusselt  number  was  6.12±0.25,  or 
an  accuracy  of  about  4  percent.  The  most  important  factor  affecting  the 
accuracy  in  this  case  was  the  accuracy  of  the  temperature  difference 
-  t^i-  The  uncertainty  calculated  here  was  the  uncertainty  based  on 
20  to  1  odds . 


Factors  Affecting  Accuracy  of  Results 

A  heat  balance  between  electrical  energy  input  and  the  thermal  en¬ 
ergy  increase  of  the  fluid  v/ill  give  some  indication  of  the  magnitude  of 
errors  (known  and  unknov/n)  v/hich  affect  reduced  quantities  such  as  local 
mixed-mean  temperature  and  v/all  heat  flux.  A  plot  of  net  electrical  en¬ 
ergy  input  (after  heat  loss  corrections)  against  measured  thermal  energy 
rise  is  given  in  figure  20.  A  deviation  from  the  line  indicates  an  im¬ 
perfect  heat  balance.  Except  for  the  two  low'est  power  runs,  the  heat 
balance  is  within  6  percent,  and  for  most  runs  it  is  much  better.  The 
heat  balance  improves  with  increase  in  povrer  and  in  almost  all  cases  the 
electrical  input  (after  loss  corrections)  is  greater  than  the  thermal 
energy  rise.  Since,  by  the  use  of  the  a-c  potentiometer,  the  electrical 
pov/er  measurement  accuracy  should  be  nearly  independent  of  power  level, 
it  is  felt  that  this  discrepancy  is  largely  due  to  errors  in  bus  bar  heat 
loss  corrections  and  to  eivors  In  measuring  bulk  te.mperature  rise.  One 
would  expect  the  percentage  error  in  measuring  a  temperature  difference 
with  a  calibrated  thermocouple  to  become  less  as  the  temperature  differ¬ 
ence  is  increased.  This  is  consistent  with  the  trend  in  the  heat  balance 
data.  It  is  concluded  that  the  unknown  error  in  overall  heat  balance  is 
small. 

There  is  some  evidence  that  there  are  unknown  errors  in  the  measure¬ 
ment  of  wall  temperature.  These  are  felt  to  be  largely  caused  by  the  in¬ 
fluence  of  stray  alternating  current  on  the  d-c  thermocouple  circuitry. 
The  thermocouples  were  spot-welded  directly  to  the  test  section,  v/hich 
was  a-c  resistance  heated.  An  electi’onic  null  indicator  was  used  with 
a  Rubicon,  Model  B,  potentiometer  to  measure  thermocouple  voltages.  It 
was  first  observed  that  the  null  indicator  lost  sensitivity  whenever  a 
test-section  thermocouple  v.'as  being  read.  When  connected  to  an  oscillo¬ 
scope,  a  60-cycle  a-c  signal  was  observed.  This  could  be  eliminated  by 
using  a  iOO-microfarad  capacitor  aci'oss  the  input  terminals  to  the  poten¬ 
tiometer  and  by  grounding  one  end  of  the  test  section  (which  previously 
had  been  floating).  These  changes  gave  good  sensitivity  except  when 
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reading  thermocouples  near  the  ungrounded  end  of  the  test  section.  This 
problem  was  corrected  by  using  a  ground  at  each  end  and  a  toggle  switch 
to  select  the  end  to  be  grounded.  The  sensitivity  was  good  at  each  end 
and  tolerable  at  the  center  of  the  test  section.  The  central  test- 
section  thermocouple  was  read  twice,  once  with  each  end  grounded.  The 
maximum  discrepancy  noted  was  0.017  millivolt  (about  0.5°  F),  and  in 
most  cases  it  was  less  than  this . 

The  foregoing  changes  are  believed  to  have  removed  the  greatest  part 
of  this  error  but  probably  not  all  of  it.  One  anomaly  was  that  the 
standardization  point  of  the  potentiometer  kept  changing  when  a  tempera¬ 
ture  traverse  of  the  test  section  was  being  made.  Apparently  the  alter¬ 
nating  current  still  getting  into  the  potentiometer  affected  the  stand¬ 
ardization  point.  The  potentiometer  was  standardized  before  each  ther¬ 
mocouple  reading. 
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forced  and  free  convection 


Sump  tank 


Mixing  box 
Approach  section 


Figure  5.  -  Assembly  photograph  of  test  and  approach 
sections  (radiation  shields  removed) . 
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Figure  8.  -  Thermocouple  circuit. 


Figure  9.  -  Comparison  of  fully  developed  Nusselt  numbers  with  analysis  of  reference 
for  upflow  with  heating. 
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Figure  10.  -  Effect  of  Rayleigh  niimher  on  entrance  Nusselt  numher  for  upflow  with  heatinj 
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Modified  dimensionless  axial  distance  from  start  of  heated  tube,  (2x/D)/RePr 

Figure  11.  -  Entrance  Nusselt  number  for  high  Rayleigh  number  run  for  upflow  with  heating. 
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Absolute  Rayleigh  number,  |Ra|,  |  (p^)3g^CpD'^A)/l6fikj 

Figure  16.  -  Comparison  of  fxilly  developed  Nusselt  numbers  with  analysis 
for  downflow  with  heating. 
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Figure  17.  -  Thermal  entrance  Nusselt  number  for  negative  Rayleigh  number  for  downflow 
with  heating. 


n  start  of  heated  tube,  x/D 
( run  ( DS ) ) . 


L  temperature  profile  for 
t  number,  1130. 


Temperal  .’.re 
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(b)  Power,  61  watts  (run  (D7)). 


Figure  IP.  -  Concluded.  Representative  wall  temperature  profile  for  down¬ 
flow  with  heating.  Peclet  number,  1130. 
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